Stability of antiphase line defects in nanometer-sized boron-nitride cones 
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We investigate tlie stability of boron nitride conical sheets of nanometer size, using first-principles 
calculations. Our results indicate that cones with an antiphase boundary (a line defect that contains 
either B-B or N-N bonds) can be more stable than those without one. We also find that doping 
the antiphase boundaries with carbon can enhance their stability, leading also to the appearance 
of localized states in the bandgap. Among the structures we considered, the one with the smallest 
formation energy is a cone with a carbon-modified antiphase boundary that presents a spin splitting 
of ~0.5 eV at the Fermi level. 

PACS numbers: 71.20.Tx, 71.15.Mb, 71.24.+q 



I. INTRODUCTION 

Boron nitride (BN) and carbon can be found in similar 
structures, from the diamond and graphite bulk forms, 
to nanostructures such as nanotub 68^*2^4 fullerenes^ 
and conical sheets iSi However, unlike carbon, BN struc- 
tures can contain three types of covalent bonds, namely 
B-N, N-N, and B-BAiS This leads to a qualitative dif- 
ference between BN and carbon in the formation of the 
graphite-derived curved surfaces that are found in nanos- 
tructures such as fuUerenes and nanotube caps. In car- 
bon, energy minimization leads to structures made of 
60° disclinations (conical sheets)ii with isolated pen- 
tagons at the apex. Carbon fuUerenes are built from such 
structures. In BN, fuUerene-like structures are formed 
from 120° disclinations. ^ These 120° BN dischnations 
are usually modeled as having a four-membered ring at 
the apexjiSiiSii^ but theoretical calculations indicate that 
non-stoichiometric BN fuUerenes, consisting of pentagon 
pairs replacing the four-membered rings, may be energet- 
ically favored, depending on the stoichiometric condition 
of growthili^ 

A possible reason for the existence of 120° disclina- 
tions in BN nanostructures, instead of the 60° ones 
that are found in carbon, is the existence of antiphase 
boundaries (APB) in BN disclinations with an odd mul- 
tiple of 60°iiS4ii The common perception is that these 
APB's, which necessarily contain a line defect of non-BN 
bonds, should present an energy cost that surpasses the 
lower elastic-energy cost of a disclination of smaller an- 
gle. However, such a "no-wrong-bond" rulei^ has been 
challenged by at least two experiments. In their work. 
Bourgeois et al& have shown the existence of BN conical 
nanostructures with a disclination angle of 300°, demon- 
strating the existence of APB's. Another example in the 
literature is the observation of disclinations in BN nan- 
otube junctions, — in experiments whose interpretation 
is also compatible with the existence of an APB. 

In the present work, we investigate the relative sta- 
bility of nanometer-sized BN cones with and without 



APB's, for N-rich and B-rich environments. Our first- 
principles results indicate that, for these two limiting 
stoichiometric situations, structures with an APB can be 
more stable than APB-free structures, consisting only of 
B-N bonds. We also find structures involving full carbon 
incorporation at the APB to be among the most stable, 
under both environments. 

The paper is organized as follows: in Sec. ^ we discuss 
the ah initio methodology for the computation of the to- 
tal energies and the definition of chemical potentials for 
the computation of the zero-Kelvin grand potentials; in 
Sec. mil we discuss our model for the formation energy 
as a sum of antiphase boundary and elastic energies; in 
Sec. llVl we discuss our first principles results for the rel- 
ative stability of the various cones and some features of 
the electronic states of the most stable structures; and 
finally, in Sec. we present some discussion of scaling to 
larger structures than considered in this work and con- 
clusions. 



II. METHODOLOGY 



A. Total energy calculations 



Our calculations are based on the density functional 
theorjilS as implemented in the SIESTA programiSS 
We make use of norm-conserving TrouUier-Martins 
pseudopotentials-^ in the Kleinman-Bylander factorized 
form,^^ and a double-C basis set composed of numerical 
atomic orbitals of finite range. Polarization orbitals are 
included for nitrogen, boron and carbon atoms, and we 
use the generalized gradient approximation (GGA)^'^ for 
the exchange-correlation potential. All the geometries 
are fully relaxed, with residual forces smaller than 0.1 
eV/A. 
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B. Zero-Kelvin grand potentials 

The comparison of BN clusters with different num- 
ber of atoms is a difficult problem. Here, we propose 
a zero-temperature thermodynamic approach based on 
the prior determination of chemical potentials to address 
this issue. The robustness of the approach will be estab- 
lished by comparison with results obtained with different 
chemical potentials, which represent different experimen- 
tal conditions, and with similar calculations in the liter- 
ature,i^.S4M*21 

In order to address the energetics of stoichiometric and 
nonstoichiometric cones, we must distinguish the differ- 
ent techniques that can be employed in their synthesis. 
In each case, we introduce suitable theoretical chemical 
potentials for nitrogen (/Xjy) and boron (/x^). For the 
case in which the cones are obtained by laser ablation or 
arc-discharge from hexagonal BN, we can have either a 
nitrogen- or a boron-rich environment, depending on the 
specific atomic reservoir employed. In the N-rich envi- 
ronment, /i^ is obtained from nitrogen in the gas phase, 
while a metallic a-B phase is used as the reservoir for the 
B-rich environment. In both cases, /ijy //^ are linked 
by the thermodynamic constraint 

I layer /-, \ 

where /i^^^"^ is the chemical potential per BN pair in the 
infinite, planar BN sheet. This constraint corresponds 
to situations in which the nanostructures are formed in 
equilibrium with BN layers. Since we are dealing with 
finite clusters, we use hydrogen atoms to saturate the 
dangling bonds at the edges. To take these additional 
H-B and H-N bonds into account, we introduce the re- 
spective chemical potentials, fJ^fjg and iJ-i^pj, and write 
the formation energy of the cones as 

Eform = Etot ~ ^bI-'-B ~ ^nI-^N ~ "^HN t^HN 

—rifjg Hug , (2) 

where Etot is the calculated total energy of the cluster, 
rig and are the number of B and N atoms, and n^j^ 
and the number of H-B and H-N bonds, re- 

spectively. The formation energy, as defined in Eq. |21 
represents the zero-Kelvin grand potential which is the 
proper thermodynamical potential for the comparison of 
the relative stability of structures with different number 
of atoms. 

A first constraint on the hydrogen chemical potentials 
is imposed by using a finite planar sheet of boron nitride 
as reference, and ascribing a null value to its formation 
energy. This allows us to write its total energy as 

E^'^'^^ — Tigpf fig-^ + rifjg Hug + n^jY ^^HN : (3) 

where n^j^ is the number of BN pairs. 

The planar reference sheet is centered at a sixfold ring 
and has Cs^ symmetry. Therefore, rigg = n^i^ and the 



preceding equation can be rewritten as 

£;f-*=nBArMBr + n«MH, (4) 

where we have defined Hh — {f^BH + Math) /2- Using 
the total energy calculations for the infinite layer and 
the finite sheet, we find /i^ = —15.46 eV. 

From the definition of fJ■H^ obtain fJ-HB it^HN) once 
f^HN (/^ffs) known. The values of one of these param- 
eters can be determined by choosing a convenient reser- 
voir. In our calculations we use the ammonia molecule 
(NH3) as a reservoir, obtaining iJ.g]\[ = —16.16 eV and 
Hhb ~ —14.76 eV in the N-rich environment, and Hhn ^ 
— 15.23 cV and Hhb = —15.69 eV in the B-rich environ- 
ment. The consistency of this methodology was tested by 
determining h^b using the BH3 molecule as a reservoir, 
and then calculating Hhn from Hh- The calculated for- 
mation energies in this case are within 0.01 eV/atom of 
those from the NH3 reservoir. For the structures contain- 
ing carbon, the chemical potentials for the carbon atom 
(Hq) and the carbon- hydrogen bonds {Hch) obtained 
from total energy calculations for bulk graphene and for 
a finite carbon sheet saturated by hydrogen atoms, re- 
spectively. 

III. FORMATION ENERGY: 
ANTIPHASE-BOUNDARY AND ELASTIC 
CONTRIBUTIONS 

The physical mechanism that defines whether struc- 
tures having APB's may be as stable as their APB-free 
counterparts is the competition between the elastic en- 
ergy associated with the cone formation and the energetic 
cost of forming an APB. Hence, we write the cone forma- 
tion energies, calculated with the above procedure, as the 
sum of two contributions. The first is the elastic energy 
required to fold a finite sheet into a cone, which we de- 
scribe by a continuum model. The second represents the 
APB chemical energetic cost associated with the non-BN 
bonds. 

The basic assumption of the continuum model is that 
the elastic energy stored in a differential element of area 
at the surface of the cone is proportional to the square of 
the surface mean curvature in that region. Within this 
model, the elastic energy of the cone is written 

= K f In + Et.piRe) , (5) 

\sina J \RcJ 

where 2a is the apex angle, R is the slant height, Rc is a 
cutoff length introduced to remove the singularity at the 
tip (linear elasticity breaks down in the neighborhood of 
the tip), and Etip{Rc) is the associated tip energy. The 
expression 

sin(.)^l-A.i_^ (6) 

relates the apex angle to the disclination angle P = N x 
60°. 
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The validity of the continuum model was first tested 
for the case of carbon, for which there is no APB. We 
found that the first-principles total energy calculations 
for the graphene sheet, and for the cones with 60° and 
120° disclinations are fitted remarkably well by a straight 
line, when plotted as a function of cos^a/ sin a, as shown 
in Fig. 1. This fact also shows that there is no tip con- 



180° and 300° (each containing an APB) were also ad- 
dressed. To make contact with experimental results for 
the synthesis of BN cones from carbon templates, we also 
considered the possibility of carbon incorporation in BN 
cones, leading to the formation of carbon zig-zag [shown 
in Fig. 2(e)] and molecular APB's. In these structures, 
every (B or N) atom along the APB is replaced by a 
carbon atom. 
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FIG. 1: Formation energy as a function of cos^a/ sin a for 
carbon cones with 60° and 120° disclinations and for the 
graphene sheet. 

tribution to the formation energy, i.e. Eup = 0, at least 
for structures with disclination angles up to 120°. This 
can be explained by the fact that these geometries have 
the general appearance of truncated cones. Figure 1 also 
allows us to consider that Rc is the same for the cones 
with 60° and 120° disclinations. 
Adding the APB energies, we write 



/ cos^a^ 






Hi) 


\ sin a J 





for the formation energy of BN cones. The constant of 
proportionality K can be obtained simply by consider- 
ing that for a 120° disclination there is no APB, hence 
^APB — ^-^^ Having K and the formation energies, we 
can use the equation above to estimate the APB energies. 



IV. FIRST PRINCIPLES RESULTS 

Figure 2 shows three of the most stable structures 
found in our study. Note that each contains a 60° discli- 
nation, leading to the formation of an APB. We consid- 
ered two kinds of APB's: one contains a sequence of par- 
allel (N-N or B-B) bonds, as shown in Figs. 2(a) and 2(c), 
and the other a sequence of zig-zag (N-N or B-B) bonds. 
These are denoted molecular APB and zig-zag APB, re- 
spectivelyii^ For cones with 120° disclinations, we con- 
sidered structures with a four-membered ring at the apex 
(without non-BN bonds), as well as structures with two 
pentagons at the apex, adjacent or not. Disclinations of 




FIG. 2: (Color online) Stable boron nitride nanometer-sized 
conical structures with 60° disclinations and an APB (see 
text). Boron, carbon, and nitrogen atoms are represented by 
white, gray, and black (blue) circles, respectively. Hydrogen 
atoms saturate the dangling bonds at the edges, (a) and (b) 
show a molecular APB of boron and the electron distribution 
of its empty gap states, respectively, (c) and (d) show the 
same for a nitrogen molecular APB and its filled gap states, 
(e) and (f) show the same for a zig-zag APB of carbon and 
its filled gap states. 

Our first principles results for the formation energy 
of cones with various disclination angles, for B-rich and 
N-rich environments, are shown in Table I. For the 60°- 
disclination cones, the third and fourth columns of Ta- 
ble I show the ratio Ej^pg/E^i, where E^i = 2.96 eV, 
from the continuum model. Note that the highest APB 
energies are associated with zig-zag configurations, more 
specifically a zig-zag boron APB in the N-rich environ- 
ment, and a zig-zag nitrogen APB in the B-rich environ- 
ment. Geometric aspects may explain these results. In- 
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deed, zig-zag APB's introduce larger strains in the struc- 
tures when compared to molecular APB's, since the par- 
allel arrangement of B-B and N-N bonds in the latter 
is more easily accommodated by the BN matrix of the 
cone. Table I also shows that for some structures the 
energetic cost of creating an APB is comparable to the 
elastic energy, i.e. Ej^pg/E^^ ^ 1. For the B-rich envi- 
ronment, these are the molecular APB of boron [shown 
in Fig. 2(a)] and the zig-zag APB of carbon (with C-B 
bonds), while for the N-rich environment, the ratio is ~1 
for the molecular APB of nitrogen [shown in Fig. 2(c)] 
and the zig-zag APB of carbon (with C-N bonds) [shown 
in Fig. 2(e)]. From these observations, we expect that 
cones with these APB's will be the best candidates to 
compete energetically with APB-free cones. 

Our results for the formation energies, shown in the 
last two columns of Table I, confirm this expectation. As 
discussed in Ref. Q, in a condition of chemical equilib- 
rium, the relative stability of the structures is determined 
by the quantity Eform/npoor^ where Upoor is the number 
of atoms of the species for which the environment is poor 
[rig (n^) in the N-rich (B-rich) case]. The underlined 
numbers in Table I indicate the most stable structures. 
A general trend that can be inferred from the numbers in 
table I is that the existence of non B-N bonds does not 
necessarily increase the formation energy. The stability 
of non-stoichiometric cones depends strongly on the en- 
vironment. This conclusion is valid for isolated non B-N 
bonds as well as for the different kinds of APB's. 

Concentrating on the structures with 120° disclina- 
tions, we note that, under N-rich conditions, the energy 
of a cone with two non-adjacent pentagons (containing 
three N-N bonds) is comparable to that of a cone with a 
four-mcmbcred ring at the apex. The same is true if the 
pentagons are adjacent. The latter is a stoichiometric 
structure with only one non-BN bond, and its formation 
energy is independent of the environment. 

Concerning the APB's in geometries with 60° disclina- 
tions, we remark again that under the two limiting sto- 
ichiometric conditions we considered, 60° disclinations 
with Ej^pg/E^i ~ 1 [the molecular APB of nitrogen (in 
the N-rich case) and the molecular APB of boron (in the 
B-rich case)] can be more stable than the APB-free struc- 
tures with 120° disclinations. These are stable structures, 
despite being responsible for the existence of gap states>i& 
(Gap states are absent in the case of cones with a 120° 
disclination.) 

The electronic charge densities associated with the gap 
states for the molecular APB's of boron and nitrogen 
are shown in Figs. 1(b) and 1(d), respectively. The 
nitrogen-APB states show a substantial degree of delo- 
calization, spreading over the BN matrix of the cone, 
while the boron- APB states remain strongly localized on 
the B atoms at the APB. In a previous study^^ we ex- 
amined the electronic levels of these two structures. In 
the nitrogen case, we found that three fully occupied lev- 
els appear in the lower half of the gap, and Fig. 1(d) 
shows that they have the same nitrogen lone-pair char- 



TABLE I: Formation energies (in eV) of stoichiometric and 
non-stoichiometric cones, for N-rich and B-rich environments. 
The first column indicates the disclination angle for each 
structure. Second column indicates either the type of APB 
or the topological defect at the apex, in each case. For the 
120° and 180° disclinations, we considered different possibil- 
ities for the defect at the apex (indicated by the number of 
squares and pentagons in each case). The third and fourth 
columns show the ratio between the the antiphase-boundary 
{Eapb) and elastic energies (-Eei) for the 60° disclinations. 
The last two columns show the formation energy per Upoor, 
as discussed in the text. Underlined numbers indicate the 
most stable structures for each environment. 
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acter that dominates the states at the top of the valence 
band of the cone. This leads to a strong mixing between 
the APB states and the valence-band states in the bulk of 
the cone, this being the reason for the delocalized charac- 
ter of the nitrogen-APB states. In the case of boron, the 
states shown in Fig. 1(b) are empty levels near the bot- 
tom of the conduction band. They have a more metallic 
character, differing from the states at the bottom of the 
conduction band. Mixing is absent in this case, leading 
to the localized distribution seen in Fig. 1(b). 

Table I also shows that cones with zig-zag APB's have 
higher formation energies, and the same occurs for struc- 
tures with large disclination angles. The result is ex- 
pected in the case of the 180° and 300° disclinations, 
since besides the cost of forming an APB, the elastic en- 
ergy is also large. However, we note that these metastable 
structures have already been found experimentallyi^ii^ 

Regarding carbon-doped cones, we find that incorpo- 
ration of carbon at the APB leads to highly stable struc- 
tures. In this case, for both kinds of APB's the most sta- 
ble structures are those where the APB is connected to 
the BN matrix by C-N (C-B) bonds in the N-rich (B-rich) 
environment. Actually, the most stable of all structures 
in our study is the C-N connected carbon zig-zag APB, 
in the N-rich environment, shown in Fig. 1(e). The elec- 
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tronic charge density of the gap states associated with 
this structure is shown in Fig. 1(f). The states are pre- 
dominantly locaHzed in the antiphase boundary, strongly 
resembling the states characteristic of the graphite tt 
band. This wavefunction localization leads also to a 
considerable spin splitting at the Fermi level. For the 
neutral charge state, the energy difference between the 
highest-occupied molecular orbital (which has majority 
spin) and the lowest-unoccupied one (which has minor- 
ity spin) is ~ 0.5 eV. This spin splitting appears only 
for structures with an odd number of electrons, a result 
that raises the interesting question about what may hap- 
pen in the limit of cones with larger sizes. Furthermore, 
we have also identified a tendency for carbon atoms to 
segregate at the APB. In our calculations, a single pair 
of carbons atoms is found to be more stable (by ^0.07- 
0.15 eV/atom) on the APB than in the bulk of a EN 
cone. This may indicate a tendency for the formation 
of these carbon- modified APB's, during the synthesis of 
BN conical nanostructures by substitution on a carbon 
template. 

Boron nitride cones can also be obtained by a substi- 
tution reactioni^ on a carbon template in which B2O3 
and gaseous N2 are used as sources of boron and nitro- 
gen, respectively. To check the stability trends in this 
case, we consider that BN cones and CO2 are the by- 
products of this reaction. The nitrogen and boron chem- 
ical potentials are then obtained in terms of the total 
energies of the species involvedn^S With this procedure, 
the value we obtained for the chemical potential differ- 
ence, — fJ'B — —189.936 eV, is very close to the value 
of -190.126 eV, obtained for an N-rich environment under 
the constraint in Eq. ^ Therefore, the stability trends 
shown in table I still hold in this case. We also checked 
the robustness of our results with respect to changes in 
the boundary conditions, which could happen, for in- 
stance, if the cones were interacting with different sur- 
faces. This was done by imposing a constant shift of 
± 0.5 eV on all hydrogen chemical potentials. The or- 
dering of the quantities Eform/npoor m table I remains 
unaltered, with relative shifts of less than 0.02 eV. Still to 
confirm the validity of our methodology, we mention that 
similar approaches have been employed in the literature 
in the determination of formation energies of defects in 
cubic boron nitride in boron nitride nanotube^ 
and in the investigation of possible structures of boron ni- 
tride fullerenesii^ In the latter work, the stability trends 



obtained by using the solid phases of boron and nitrogen 
together with the constraint in Eq. ^ remained largely 
unchanged when the authors considered the gas phase 
limit of atomic reservoirs. 



V. REMARKS ON SIZE SCALING AND 
CONCLUSION 

Finally, we discuss the relevance of our results in face 
of the expected scaling trends. As the size of a cone 
increases, the APB energy scales linearly with the slant 
height R (see Eq. O, i.e. with the square root of the 
total number of atoms in the structure, while the elastic 
energy scales sub-linearly with R, as shown in Eq. \^ 
Hence, structures without APB's become energetically 
favored for sufficiently larger sizes. However, cones of 
the size we consider in this work may represent the seeds 
present in the initial steps of the nucleation process that 
generates larger cones or other nanostructures. Not only 
the energetics, but also kinetics play an important role in 
this process, and to the extent that we may expect low 
energy configurations to be favored for the seeds, cones 
with certain types of APB's may be good candidates, 
and kinetics may determine the formation of larger size 
structures containing APB's. 

To summarize, we used first-principles calculations to 
address the question of how stable are antiphase bound- 
aries in boron nitride nanometer-sized conical nanostruc- 
tures. Our results indicate that 60° disclinations with an 
APB can be more stable than APB-frec 120° disclina- 
tions. We show that the formation energies of the stable 
APB's are comparable to the disclination elastic energy, 
and that the dependence of the latter on the disclina- 
tion angle is very well described by a continuum model. 
Also, we observe that carbon incorporation at the APB 
leads to highly stable structures, of which the most sta- 
ble presents a large spin splitting (~ 0.5 eV) at the Fermi 
level. 



Acknowledgments 

We acknowledge support from the Brazilian agen- 
cies CNPq, FAPEMIG, and Instituto do Milenio em 
Nanociencias-MCT. 



* Electronic address: |chacham@fi sica.uftng . jpr | 
^ S. lijima. Nature (London) 354, 56 (1991)7 
^ N. G. Chopra, R. J. Luyken, K. Cherrey, V. H. Crespi, 

M. L. Cohen, S. G. Louie, and A. Zettl, Science 260, 966 

(1995). 

^ Z. Weng-Sieg, K. Cherrey, N. G. Chopra, X. Blase, Y. 
Miyamoto, A. Rubio, M. L. Cohen, S. G. Louie, A. Zettl 
and R. Gronsky, Phys. Rev. B 51, 11229 (1995). 



* A. Loiseau, F. Willaime, N. Demoncy, G. Hug, and H. 

Pascard Phys. Rev. Lett. 76. 4737 (1996). 
5 H. W. Kroto, J. R. Heath, S. C. O'Brien, R. F. Curl, and 

R. E. Smalley, Nature (London) 318, 162 (1985). 
® D. Golberg, Y. Bando, K. Kurashima, and T. Sasaki Appl. 

Phys. Lett. 72, 2108 (1998). 

A. Krishnan, E. Dujardin, M. M. J. Treacy, J. Hugdahl, S. 
Lynum, and T. W. Ebbesen. Nature 388 451 (1997) 



6 



L. Bourgeois, Y. Bando, W. Q. Han and T. Sato. Phys. 

Rev. BGl, 7686 (2000). 
^ P. Zhang and V. H. Crespi. Phys. Rev. B 62 11050 (2000) 
° T. Dumitrica, H. F. Hettinger, G. Scuseria and B. I. Yakob- 

son. Phys. Rev. B 68 085412 (2003) 
^ For a definition and examples of disclinations in BN nanos- 

tructures see referencei& 
^ Simone S. Alexandre, Mario S. C. Mazzoni, and Helio 

Chacham Appl. Phys. Lett. 75, 61 (1999). 
^ F. Jensen and H. Toftlund. Chem. Phys. Lett. 201, 89 

(1993). 

* Simone S. Alexandre, Helio Chacham, R. W. Nunes. Phys. 
Rev. B 63 045402 (2001). 

^ P. W. Fowler, K. M. Rogers, G. Seifert, M. Terrones, and 
H. Terrones, Chem. Phys. Lett. 299, 359 (1999). 

^ S. Azevedo, Mario S. C. Mazzoni, H. Chacham, R. W. 
Nunes. Appl. Phys. Lett. 82 2323 (2003) 

^ A. Loiseau, F. Willaime, N. Demoncy, N. Schromchenko, 
G. Hug, C. CoUiex and H. Pascard Carbon. 36. 743 (1998). 

* D. Golberg, Y. Bando, L. Bourgeois, K. Kurashirna, and 
T. Sato, Appl. Phys. Lett 77, 1979 (2000). 

^ W. Kohn, L. J. Sham. Phys. Rev. 140, A1133 (1965). 

'° D. Sanchez-Portal, P. Ordejon, E. Artacho and J. M. Soler. 

Int. Journ. of Quant. Chem. 65, 453 (1997). 

N. Troullier and J. L. Martins. Phys. Rev. B. 43, 1993 

(1991). 

L. Kleinman and D. M. Bylander. Phys. Rev. Lett. 48, 



1425 (1982). 

J. P. Perdew, K. Burke, M. Ernzerhof. Phys. Rev. Lett. 77 
3865 (1996) 

W. Orellana and H. Chacham. Appl. Phys. Lett. 74, 2984 
(1999) 

W. Orellana and H. Chacham. Phys. Rev. B. 62, 10135 
(2000) 

W. Orellana and H. Chacham. Phys. Rev. B. 63, 125205 
(2001) 

T. M. Schmidt, R. J. Baierle, P. Piquini, and A. Fazzio. 
Phys. Rev. B. 67, 113407 (2003). 

There is a very small diflerence in slant height between 
our 60°- and 120°-disclination structures, which leads to a 
negligible change in the logarithmic term in Eq. 15] between 
these two types of structures. This has been neglected in 
our determination of K. 

The substitution reaction described allows us to write the 
nitrogen and boron chemical potentials as fj,N = \E{N2) 
and /is = \[E{B203) + '^E{C) - |£;(C02)1, where E(N2) 
and E(C02) are the total energies of the nitrogen and car- 
bon dioxide molecules, respectively, E{B203) is the total 
energy of a planar hexagonal B2O3 structure and E(C) is 
the energy per atom in graphene. All these quantities are 
obtained by means of first-principles calcultions and the 
resulting values of fiN and /is are -270.151 eV and -80.214 
eV, respectively. 



